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Abstract — This article presents a detailed overview and as- 
sessment of contact-free UV light discharge systems (UVDS) 
needed to control the variable electric charge level of free- 
flying test masses which are part of high precision inertial 
sensors in space. A comprehensive numerical analysis approach 
on the basis of experimental data is detailed. This includes 
UV light ray tracing, the computation of time variant electric 
fields inside the complex inertial sensor geometry, and the 
simulation of individual photo-electron trajectories. Subsequent 
data analysis allows to determine key parameters to set up 
an analytical discharge model. Such a model is an essential 
system engineering tool needed for requirement breakdown and 
subsystem specification, performance budgeting, on-board charge 
control software development, and instrument modeling within 
spacecraft end-to-end performance simulators. Different types 
of UVDS design concepts are presented and assessed regarding 
their robustness and performance. Critical hardware aspects like 
electron emission from air-contaminated surfaces, interfaces with 
other subsystems, and spacecraft operations are considered. The 
focus is on the modeling and performance evaluation of the 
existing UVDS on board LISA Pathfinder, an ESA technology 
demonstrator spacecraft to be launched in 2014. The results 
have motivated the design of a more robust discharge system 
concept for cubical test mass inertial sensors for future space 
missions. The developed analysis tools have been used for 
design optimization and performance assessment of the proposed 
design. A significant improvement of relevant robustness and 
performance figures has been achieved. 
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Nomenclature 

AC = Alternating current 

AIT = Assembly, integration, and test 

BRDF = Bidirectional reflectance distribution function 

CAD = Computer-aided design 

DC = Direct current 

DFACS = Drag-free and attitude control system 

DLR = German aerospace center 

DoF = Degree of freedom 

EH = Electrode housing with electrodes 

ESA = European Space Agency 

GCR = Galactic cosmic ray 

LED = Light-emitting diode 

LISA = Laser Interferometer Space Antenna 

EPF = LISA Pathfinder 

LTP = LISA Technology Package 

MEW = Measurement bandwidth 

NGO = New Gravitational Wave Observatory 

PDF = Probability density function 

PTE = National Metrology Institute of Germany 

SEP = Solar energetic particle 

SC = Spacecraft 

TM = Test mass 

UV = Ultra-violet 

UVDS = UV light discharge system 

UVLF = UV light feedthrough 

I. Introduction 
A. Problem Formulation 

INERTIAL sensors are key components of scientific space 
missions aimed at measuring the effects of spacetime cur- 
vature caused by celestial bodies. Spacetime curvature can be 
measured directly using at least two free-falling test masses as 
gravity references in space, by observing the distance variation 
between these reference points. The variation between a test 
mass and its hosting spacecraft is used as reference for the 
spacecraft drag-free control system ||2l|[3l needed to suppress 
the spurious non-gravitational forces acting on the shielding 
spacecraft. 

The envisaged space-born gravitational wave detector 
LISA/NGC|^|4 ||5 | relies on the direct measurement of space- 
time curvature, using laser interferometry for the distance 
measurement and high precision inertial sensors with cubical 

^Free-falling particles are not disturbed by any forces and move along their 
geodesic lines in curved space-time 1 1 1. 

^Since 2011, in the framework of an ESA only mission, LISA is called 
NGO. 
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test masses as the gravity references. A description of a 
modern high precision inertial sensor for space is given in 
Appendix |A| The currently assembled LISA Pathfinder (LPF) 
spacecraft is an ESA/NASA technology demonstrator (to be 
launched in 2014 and operated in an orbit around the LI 
Lagrange point), aimed at testing novel LISA technologies 
161 . Its scientific pay load, the LISA Technology Package Q 
is, among other instruments, equipped with two high precision 
inertial sensors. The sensors must demonstrate a level of 
residual, non-gravitational acceleration noise in the order of 
10~^^ m/s^/Hz^/^ or below, down to millihertz frequencie^ 

The cubical test masses of the LPF and the LISA/NGO 
inertial sensors have no mechanical contact with their sur- 
rounding housing structures. Ideally, in the absence of any 
non-gravitational forces, the test masses should move along 
geodesic lines. However, in the real instruments, pure geodesic 
test mass motion is impaired by non-gravitational forces |8| 
that accelerate the test masses away from their geodesic lines. 

A major source of non-gravitational disturbance noise is 
caused by the net electric charge Q^^ accumulating on the 
test masses. Sources of test mass charging are: 

1) Contact electrification O due to test mass release 
into free-flight from the release mechanism (see Ap- 
pendix |A]). 

2) Energetic, charged particles such as galactic cosmic rays 
(GCRs) and large amounts of high-energy plasma, peri- 
odically expelled by the sun in solar energetic particle 
(SEP) events. These particles penetrate the spacecraft 
shielding and eject or deposit charges on the test mass 
|10||11|. 

Taking LISA Pathfinder as an example, from the first effect, 
it is expected that a charge level up to ±2.35 • 10^ elementary 
charges (e) with arbitrary sign will deposit on a test mass 
|[T2ll . This corresponds to a test mass potential Vtm = ±1.1 V 
relative to the grounded housing surfaces. From the second 
type of effects, a charge rate of about ±50 e/s is expected due 
to GCRs for solar minimum conditions which corresponds to 
approximately ±4.3-10^ e (±20 mV) within 24 h. Furthermore, 
about ±5 • 10^ elementary charges accumulate due to one 
typical, small SEP event 1 13 |. The probability of such an event 
occurring in a 24 h period is estimated to be 2.4 % (calculated 
from the data reported in (131 Fig- 3] for the year 2011). 

Thus, 24 h after test mass release, up to ±2.39 • 10^ e 
(Vtm = ±1.12 V) may deposit on the test mass — assuming a 
positive charge level has been accumulated on the test mass 
due to contact electrification after release. With the same 
argument, after six month without re-grabbing the test mass, 
a potential of up to ±4.9 V might be present on the test mass 
(assuming 4 smaller SEP events). 

Figure [T] shows the contribution of all acceleration noise 
effects on the test mass at 1 mHz, caused by a coupling 
with the test mass charge Qtm- The linear spectral density 
Sa,QTM of the acceleration noise at ImHz is plotted versus 

^The LISA Pathfinder requirement on differential acceleration noise be- 
tween two test masses along the sensitive axis is 3-10"-'^'^ m/s^/Hz-*^/^ in the 
measurement bandwidth (MBW) from 1 — 30 mHz. The LISA acceleration 
noise requirement is 3 • 10~^^ m/s^/Hz^/^ in the MBW from 0.1 - 100 mHz. 
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Fig. 1. Linear spectral density of acceleration noise effects due to coupling 
with accumulated test mass charge at ImHz. The budget is shown for the 
reference inertial sensor described in Appendix [A| 

increasing test mass potential Vtm = Qju/^tou where Ctot is 
the total electrostatic capacitance between the test mass and 
its surrounding housing. Sa,QTM is the root- sum- square value 
of multiple contributors that couple with Qj^ (e.g., coupling 
of voltage noise from electrostatic actuation, stray voltage 
fluctuations, increase of stiffness coupling with the control 
jitter, etc.Q Note that the linear spectral density of the total 
residual test mass acceleration noise is given by the sum of the 
charge related effects and various other contributors |8||14|. 
However, Figure [T] shows that with increasing test mass charge, 
also Sa,QTM increases, such that, at Vtm = 900 mV, the 
acceleration noise at ImHz is already at 10~^^ m/s^/Hz^/^ — 
only due to test mass charge related effects. 

Since the LISA Pathfinder mission goal is to properly 
characterize all effects causing acceleration noise below 
lO-^^^m/s^/Hz^/^the test mass charge has to be much smaller 
to limit its acceleration noise contribution. Moreover, the 
LISA/NGO requirement on the residual test mass potential 
is 0.5 mV (Qtm = 1 * 10^ ^) order not to disturb the science 
measurements. The verification of the required LISA/NGO 
charge level control accuracy is one of the experiments to be 
demonstrated with LISA Pathfinder ifTSl . 

Thus, in order to meet the stringent acceleration noise 
requirements, positive and negative charges must be removed 
from the inertial sensor test masses ("bipolar discharging"). 
Charge control by using a thin conducting gold wire as an 
electrical connection between the test mass and the electrode 
housing structure (e.g., applied in modem accelerometers 
as used in GOCE |[16]| or MICROSCOPE (TT\) is not an 
option for inertial sensors with residual acceleration noise 
requirements in the order of 10~^^ m/s^/Hz^/^. According to 
1 18 1, such wires introduce increased stiffness and damping. 
The damping force itself is negligible, but the thermal noise 
it causes (fluctuation-dissipation theorem) limits the perfor- 

Additional acceleration noise effects at ImHz due to charged particle 
impacts from the environment and charge control actuation are not considered 
in Sa,QTM since it contains only effects that couple with accumulated DC test 
mass charge Qtm- 
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mance of an LPF and LISA/NGO like inertial sensor to 
4.5 • 10"^^m/s^/Hz^/^ at ImHz (i.e., already above the total 
noise requirement). 

As a consequence, the test mass charge has to be controlled 
without any mechanical contact to the surrounding housing 
structure. This can be achieved by using ultra-violet (UV) light 
via the photoelectric effecj^ 

B. Previous work and significance this article 

A contact-free, photoemission based UV discharge system 
has been used to control the charge of the four electrostatically 
suspended spherical test masses on-board Gravity Probe B 
120 J . The surfaces of the test mass and of a dedicated charge 
control electrode (integrated in the surrounding housing struc- 
ture) have been illuminated with UV light; the direction of the 
induced photo-currents from both surfaces has been controlled 
by application of either positive or negative DC bias voltages 
to the charge control electrode (3V with respect to the test 
mass surface). 

An early design of the discharge system for the cubical, gold 
coated test masses of the LISA Pathfinder and LISA inertial 
sensors is described in [21 1. A modification of the UV light 
injection angle into the sensor from the original design is de- 
scribed in |22|. An analysis of the light injection modification 
and the influence on achievable discharge performance and 
robustness is not reported therein. 

In |23| a simulation of closed-loop drag-free and attitude 
control including charge control operations is mentioned; the 
description of the used discharge system model is omitted. In 
(241 a simplified discharge model is reported, similar to the 
model previously developed by the authors [25 1. The simpli- 
fied model treats specific pairs of housing surfaces and their 
adjacent surfaces on the test mass as parallel plate capacitors, 
where a uniform electric field (neglecting edge effects) is 
assumed. Photo-electron emission is assumed perpendicular 
to the emitting surface and therefore constrains the released 
electrons to the region between the two adjacent surfaces. 
However, a major part of the electron emission within the 
LPF inertial sensor happens at comer regions, where the 
illuminated housing parts have no adjacent test mass surfaces 
and the electrical fields can not be described by parallel plate 
capacitors. 

Discharge concepts using UV LEDs are described in 
various publications (e.g., 1261127 1). However, a design 
concept for the cubical, gold coated test-masses of the LPF 
and LISA/NGO inertial sensors as well as a model-based 
design justification including robustness and performance 
assessments (e.g. considering the effects of photoemission 
from air-contaminated surfaces, UV light ray tracing in 
complex sensor geometries, and calculation of electron 
transition ratios between relevant surfaces) is not presented. 

This article summarizes the operational principle of UV 
light discharge systems and gives an overview of different 

^Also alternative techniques for the implementation of contact-free test 
mass discharge systems have been discussed. These techniques make use of 
radioactive sources, ion sources, or field emission cathodes (see e.g. 1191). 



design concepts in Section [ll| A comprehensive toolbox for the 
mathematical modeling of such discharge systems is presented 



in Section III It considers crucial aspects like photoemission 
from air-contaminated surfaces with realistic models for the 
kinetic energy distribution and the angular distribution of the 
photo-electrons, the implementation of detailed geometrical 
sensor features, electron propagation through complex and 
time-variant electric fields that have been computed with a 
finite-element tool, as well as UV light scattering measure- 
ments from rough surfaces. The toolbox is used to derive an 
analytical model, needed for the breakdown of discharge rate 
level requirements down to requirements on subsystem level. 
Such a model is also essential for the development of the on- 
board charge control software, for the simulation of closed- 
loop charge control performance by means of spacecraft end- 
to-end simulators (see 1281 ), and for performance budgeting 
and control during the spacecraft integration and verification 



process. In Section |IV| the developed tools are used to analyze 
the LPF discharge system design ll2T1lll22ll and to justify 
important design modifications. 

Aided by the toolbox, the robustness of different design 
concepts can be conveniently evaluated. A more robust, model- 
based design for future space inertial sensors is presented 
in Section |V| and also evaluated in terms of robustness and 
performance. The design concept makes use of synchronized, 
high-frequency switching of UV LEDs, combined with an 
ideal light injection into the inertial sensor. In this way, 
the already existing high-frequency voltages (applied by the 
injection electrodes for electrostatic measurement purposes) 
will ideally assist robust and bipolar test mass discharging. 

II. Overview of UV Light Discharge Systems 

This sections describes the basic principles of contact-free 
UV light discharge systems for inertial sensors and introduces 
different design concepts. Critical aspects like the sensitivity of 
the photoemission process on air-contaminated sensor surfaces 
|29 1 and the UV wavelength are addressed. 

A. Principle of Contact-Free UV Light Discharge Systems 

The principle to discharge a free-flying test mass using UV 
light is illustrated in Figure [2] The basic idea is to illuminate 
specific surfaces of the test mass and the surrounding housing 
with UV light such that photo-electrons are released from 
these surfaces via the photoelectric effect. A positive test mass 
discharge rate Qtm > is obtained when the electron flow 
from the test mass to the housing dominates the electron flow 
from the housing to the test mass. A negative discharge rate 
Qtm < is obtained when the electron flow from the housing 
to the test mass is dominating. The UV light is injected via UV 
light feedthroughs (UVLF). Usually, a significant amount of 
light is reflected between individual inertial sensor surfaces. 
Thus, in addition to the primary illuminated surfaces, also 
other surfaces absorb UV light. The fraction of injected light 
absorbed by an arbitrary surface i (taking multiple reflections 
into account) is denoted pi . In Figure [2j the totally absorbed 
light on the test mass is indicated by the shaded area denoted 
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Fig. 2. Principle of test mass discharging shown for a negative discharge rate. 
The key parameters are the absorbed UV light on EH and TM, the amount 
of electrons released per absorbed UV light (quantum yield), and the photo- 
electron transition ratios between EH and TM surfaces. The transition ratios 
are strongly influenced by the electric fields and the initial velocity vector vq 
of the photo-electrons. 



PxestMass- Only UV light that is absorbed by a surface con- 
tributes to the photoelectric effect. 

In principle, when the energy of the absorbed light equals 
the material work function, photoemission occurs for those 
electrons occupying the highest energy level (Fermi). Increas- 
ing the photon energy even further will also lead to emission 
of valence electrons below the Fermi level, which then have 
sufficient energy to be excited above the vacuum level 1301 . 
According to Fowler's derivation (STJ, the photo-current / for 
zero Kelvin should scale around the emission threshold (i.e., 
when electrons are removed from energy levels close to the 
Fermi energy) as 



/ = A-ihiy-^-eVif 



(1) 



for {hu -<!>)> and {hu - ^ - eVi) > 0. In Eq. [T| 
is the photon energy, A is a constant, ^ the material work 
function, and Vi an arbitrary voltage that is applied between the 
emitting and the receiving surface. Hence, the photo-current 
from one surface to another constitutes the total of all emitted 
electrons up to a certain energy of the voltage Vi. The electric 
field caused by the voltage Vi enhances or attenuates the 
electron flow. V originates from AC and DC voltages applied 
to individual electrodes (see eqs. 1_I-13 in Appendix [A|, and 
from the test mass potential Vtm- The distribution of kinetic 
energies normal to the material surface is approximatel}}^ 
given by the derivative of Eq. [T] The exact kinetic energy 
distribution of the photo-electrons is, however, obtained from 
the analysis of dedicated measurement campaigns |29|. The 
photo-current between representative sample surfaces and a 
collecting anode has been measured while, a bias voltages 
was stepwise increased until the photo-current was entirely 
suppressed. In this way, the total energy distribution can be 
derived; it includes important features like the distribution at 
about 300 K (which is within the operating temperature range 
of the inertial sensor) and the low-energy modification due 
to energy dependent transmission barriers. A typical kinetic 

^Equation and Fowler's law apply to the normal energy only; the total 
energy distribution is given by DuBridge's equation, according to 1321 . 



energy distribution (derived for gold surfaces under 253.6 nm 
UV illumination) is shown in Section |III-B| 

Whether individual surfaces contribute to positive or neg- 
ative test mass discharge rates depends on the trajectories of 
the photo-electrons in the volume between the test mass and 
the housing structure. In addition to the present electric fields, 
the trajectories depend on the position, the direction, and the 
kinetic energy of the photo-electrons. In Figure |2j the initial 
state of the photo-electrons is indicated by the velocity vector 

The number of released electrons from a surface per ab- 
sorbed light is characterized by the quantum yield. The total 
quantum yield comprises contributions from the bulk material 
on the one hand and effects relating to the material surface 
on the other hand. The photo-current of Eq. [T] only takes 
bulk emission into account. Phenomena of the material surface 
(e.g., molecular adsorption and associated dipole moments) 
enter this expression only through their effect on the work 
function as a transmission barrier which determines the width 
of the kinetic energy distribution and consequently the yield. 
On the other hand, surface effects, such as the vectorial photo- 
effect (331, are highly sensitive to polarization and incidence 
angle of the light which enables their experimental separation 
from bulk emission. Surface effects have been found to in- 
crease the yield up to 20 % at incidence angles of around 60° 
with respect to the surface normal |29|. The total quantum 
yield of an arbitrary surface is illustrated by the number of 
emitted electrons in Figure [2] 

B. Photoemission Sensitivity 

The quantum yield of gold surfaces under UV illumination 
with energy close to the surface work function is heavily 
influenced by gross contamination from packaging, shipping, 
and handling, as well as by contaminant substances adsorbed 
from the atmosphere, in particular water and hydrocarbons. 
Although gross contamination is removed from the discharge 
relevant surfaces by an initial (plasma) cleaning step right 
before the inertial sensor assembly, integration, and test (AIT) 
procedure, it can not be avoided that these surfaces will be 
exposed to air before they are brought into vacuum with a 
typical initial pressure of 10~^mbar. 

As shown in 12911 , the emission efficiency of cleaned (Ar- 
lon sputtered) gold surfaces initially increases due to the air 
exposure, before diminishing below the detection limit of 
the used measurement apparatus on a timescale of several 
hours. Outbaking at about 125 °C, restored the emissivity 
and also reduced the quantum yield variation below a factor 
of approximately three (calculated as the ratio between the 
minimum and maximum quantum yield from 8 independent 
sample measurements). 

If the gold surfaces are contaminated, their work function 
changes to ^' = ^ — due to the adsorption of the 
contaminants, where is typically between 0.3 V and 0.8 V 
|29|. Assuming the photo-current obtained from bulk emission 
scales according to Eq. [T] this implies that a work function 
variation will change the photo-current dramatically if hu 
is very close to ^, but only moderately if hu ^.In order to 
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reduce the sensitivity to surface contaminations, it is therefore 
desirable to have hu ^ so that surface contaminations 
have only a moderate impact on the total photo-current. Work 
functions for air-exposed gold surfaces are typically in the 
order of ~ 4.1 eV 13411351 , which we also found from initial 
measurements with minimal water adsorption ||29J . 

C. UV Light Actuation Strategies 

As indicated in Figure [2j the electric fields strongly affect 
the electron trajectories and, obviously, they can be used to 
enhance or attenuate the transition of electrons. In combination 
with the selected UV light injection into the inertial sensor, 
this effect can be utilized for test mass discharging purposes. 
In this context, one may distinguish two different UV light 
actuation concepts: 

1) Switch on the UV light only at periods in time when 
the present voltages support the desired electron flow 
and consequently attenuate (or suppress) the unwanted 
electron flow. Ideally, for positive test mass discharge 
rates, the electron transition ratios from the TM to the 
EH surfaces are equal to one (i.e., all electrons make the 
transition) and those from the EH to the TM surfaces 
are equal to zero (i.e., no electron makes the transition), 
and vice versa for negative discharge rates. The concept 
works although the test mass and the electrode housing 
surfaces are illuminated at the same time (as long as 
both sides emit photo-electrons). 

2) Switch on the UV light constantly, irrespective of any 
applied voltages. Ideally, a positive discharge rate is 
obtained when only the TM is illuminated; a negative 
discharge rate for pure EH illumination. Usually, both 
sides are illuminated due to reflections; therefore, bipolar 
discharging is obtained by "adjusting" a specific imbal- 
ance of EH and TM illumination (achieved by different 
UVLFs which either point towards the EH or towards 
the TM). The transition ratios are somewhere between 
and 1, depending on the (time- variant) voltages applied 
to the electrodes. 

The first concept is referred to as "synchronized AC charge 
control" in the following, and the latter as "unsynchronized 
DC charge control". In this context, AC and DC refer to the 
on/off switching frequency of the UV light with respect to 
the applied voltages. AC switching means to turn on the light 
only for a short fraction of time (compared to the period of the 
alternating voltages applied to the electrodes). DC switching 
means that the light is turned on much longer, such that it can 
be considered constant w.r.t. the alternating voltages. Ideally, 
the second concept generates bipolar test mass discharge rates 
without application of any voltages; the averaged transition 
ratios can be increased or decreased through application of 
DC bias voltages. The first concept requires the application 
of alternating (e.g., sinusoidal) voltages to obtain bipolar 
discharge rates. 

III. Discharge Model 




Fig. 3. Visualization of the three-dimensional geometry model of the LPF 
inertial sensor as used for ray tracing simulations. The scene shows a test mass 
corner with the spherical launch lock interface (upper part of the figure) inside 
the housing structure with the electrodes and caging finger holes (bottom part). 
The scene is illuminated by the UVLF pointing towards the test mass (no. 1 
in Fig. [Is) . Also visualized are some UV light rays and the absorbed light 
power at the impact positions. 



analyze UV discharge systems. Dedicated hardware measure- 
ments can be taken into account as model input parameters. 



The following sections [III- A| and [III-B | describe comprehen- 
sive numerical tools, which have been developed to model and 



Section III-C introduces an analytical discharge model which 
makes use of outputs from the developed analysis tools. 

A. Ray Tracing Tool 

In order to predict photo-currents from individual inertial 
sensor surfaces, the absorbed UV light power of these 
surfaces must be known. The spatial light absorption within 
an integrated sensor cannot be measured directly. For simple 
geometries, the light distribution may be calculated using 
analytical methods; however, this is not possible for the 
complex inertial sensor geometry. Therefore, a ray tracing 
tool has been developed to calculate and visualize the 
absorbed light power of arbitrary sensor surfaces. The tool 
approximates the exact solution by tracing a large number of 
generated rays through the scene. When a ray hits a surface 
of the geometry, it is partly absorbed and partly reflected, 
according to the material's reflection and scattering properties. 
The reflected ray is further traced through the scene until 
its remaining power is negligible. Figure [3] shows the gap 
between test mass and electrode housing, generated by the 
three-dimensional geometry model of the reference inertial 
sensor (see Appendix |A]). Also the ray-tracing results are 
visualized by a limited number of UV light rays together with 
the absorbed light power at the corresponding light impact 
positions. 

1 ) Ray tracer inputs and measurements: The scene geom- 
etry is directly imported from the CAD model of the inertial 
sensor flight hardware. Hence, all mechanical design details 
(e.g., specific electrodes, caging finger holes, test mass cor- 
ners, etc.) are completely considered in the three-dimensional 
geometry model. The materials of all surfaces can be assigned 
since they define the physical properties as needed by the ray 
tracing simulation. 

A light source model generates the light rays as emitted 
from the tips of the UVLFs. The rays are randomly generated 
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Fig. 4. Processed measurement data of a typical two-dimensional UV light 
intensity distribution as emitted from an UVLF with a 1 mm multimode fiber 
inside. The raw data used to generate the distribution has been provided 
by Imperial College London, the manufacturer of the LPF UV light source 
hardware. 
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Fig. 5. Reflection curve of a typical electrode gold surface (Ra < 10 nm) 
versus UV light illumination with 254 nm wavelength at different incident 
angles. The measurement of the non-polarized reflection curve has been 
performed by ESA on a representative surface sample of an inertial sensor 
electrode. 



such that the simulated Hght source characteristics corresponds 
to the wavelength spectrum and the spatial intensity distribu- 
tion of the real light source. As an example, the processed 
measurement data of a typical light intensity distribution from 
an UVLF tip is shown in Figure |4] The distribution is one of 
the input parameters needed by the random number generator 
of the light source model. 

Each time a light ray hits an inertial sensor surface, some of 
its power is absorbed by the surface and the rest is reflected. 
The amount of reflected and absorbed light depends on the 
material, the angle of incidence (defined with respect to the 
surface normal), the wavelength and the polarization of the 
light. The reflectance of a typical inertial sensor gold surface 
versus angle of UV light incidence is shown in Figure [5] 

For the scattering of the reflected light, a simple rule 
of thumb is that specular reflection (light incident angle = 
reflection angle) can be assumed if the mean surface roughness 
Ra is smaller than 1/20 of the light wavelength A. Otherwise, 
diffuse scattering has to be considered to obtain a realistic 



light propagation after reflection from "rough" surfaces. The 
amount and distribution of UV light scattering from rough 
sensor surfaces is characterized by the bidirectional reflectance 
distribution function (BRDF). The BRDF is a 5 -dimensional 
function and describes the ratio of the reflected radiance 
exiting from the surface in a particular direction, to the directed 
irradiance incident on the surface. It is defined as 

BRDF(^„(/)„^„(/)„A) = J"' (2) 

where the subscripts i and s mean incident and scattered. Vts 
is the solid angle at which the scattered radiance is measured 
and depends on the size of the detector; the term cos(6>5) 
is a correction factor to adjust the illuminated area to its 
apparent size when viewed from the scatter direction. The 
angles 0^ and 6s are the azimuth and the zenith angle with 
respect to the surface normal. They describe the hemisphere 
above the surface in which the radiance P^ of the reflected 
beam can be characterized. The BRDF further depends on 
the azimuth angle of the incident beam (only needed for 
anisotropic materials), its zenith angle 9i (corresponding to 
the incidence angle w.r.t the surface normal), and on the 
light wavelength A. In the ray tracer, BRDF functions are 
represented by at set of measurements, obtained from a 
robot-based gonioreflectometeiQ apparatus 13611 for various 
light incident angles. Practically, the BRDF can only be 
measured for a limited combination of the five relevant 
parameters; therefore, a weighted linear interpolation method 
is used to derive data sets for all possible light incident 
angles and wavelengths from the available measurement sets. 
Subsequently, for each incident angle, a two dimensional 
light distribution probability density function (PDF) can 
be calculated from the interpolated measurement data. The 
direction of the reflected light ray is then randomly generated 
according to this PDF. 

2 ) Ray tracer outputs: The output of the ray tracer is a list 
of all occured ray impacts on the inertial sensor surfaces. For 
each impact, the absorbed light power, the impact position, 
the angle of incidence with respect to the surface normal, 
and the light wavelength of the incident ray are stored. The 
integrated power absorption of all impacts on an arbitrary 
surface i (e.g., an actuation/sensing electrode, the test mass, 
etc.) is denoted illumination ratio pi, where i is the surface 
identification number. The illumination ratios of all surfaces 
can be derived by post-processing of the ray-tracing data 
and correspond to a fraction of the total injected light power 
such that Pi = 1. The accuracy of the solution increases 
with the number of generated rays and converges towards the 
analytical solution for an infinite number of rays. For a typical 
inertial sensor illumination scenario, convergence of relevant 
surface illumination ratios (i.e., the numerical variation of the 
calculated values is less than 1%; looking only at surfaces 
where pi > 1.0%) is reached by tracing ^ 10^ rays. 

^In this context, a gonioreflectometer is a device which allows the precise 
control of the angles of the incident and reflected light beams in a reflection 
measurement. 
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B. Electron Tracing Tool 

In order to predict discharge rates between inertial sensor 
surfaces, the probabiHty of released photoelectrons to travel 
from surface i to another surface j must be known. The 
prediction of individual electron trajectories is sophisticated, 
especially for complex geometries (e.g., like the illuminated 
caging finger holes and the test mass corner spheres) and 
in regions with inhomogeneous electric fields. Note that the 
LISA Pathfinder discharge system design strongly illuminates 
exactly such regions (see Figure [Sj. 

For this reason, an electron tracing tool has been developed 
to propagate the trajectories of individual electrons from their 
origins to their absorption points on the receiving surface. The 
complex and time-variant electric field in the gap between 
the surfaces is considered when the electron trajectories are 
simulated. The goal is to obtain the transition ratios between 
all pairs of surfaces for a specific interval of time. 

1 ) Electron tracer inputs and measurements: The electron 
trajectories are simulated considering the following inputs: 

• CAD model of the inertial sensor geometry 

• Absorbed light power and impact position as simulated 
by the ray-tracing tool 

• Measured kinetic energy probability density function of 
the photo-electrons 

• Angular distribution of the photo-electrons 

• Time-variant electric fields due to applied voltages and 
test mass charge 

• Time of electron release w.r.t. the applied voltages 

• Measured quantum yields of the emitting surfaces 

For each electron, an initial condition (time, release position, 
and velocity vector) is needed. Based on these initial condi- 
tions and the instantaneous electric field, electrons are traced 
through the geometry by solving their equations of motion. 
The tracing is stopped when the electron hits its receiving 
surface. 

The release position of the electrons is obtained from the 
light impact positions as previously calculated by the ray 
tracing tool. 

The initial velocity is computed from the kinetic energy 
distribution of the photo-electrons. A typical distribution, used 
as electron tracer parameterization for the gold surfaces, is 
shown in Figure [6] It closely resembles the distribution for 
the total energy according to DuBridge |32| for emission 
from a free electron gas at 300 K; which has been found by 
detailed analysis of measurement data obtained from specif- 
ically conducted sample measurement campaigns 1291 . The 
analysis considers the electrode geometry of the measurement 
apparatus, the impact of dipolar adsorbants, and the effect 
of disturbing electric fields that have been observed during 
the measurements conducted at the German Aerospace Center 
(DLR). For each electron to be traced, a random initial kinetic 
energy Ekin(Xi) is generated by mapping a pseudo random 
number Xi between and 1 (generated from a uniform 
distribution) to the input range of an inverse cumulative 
distribution function. The cumulative distribution is derived 
from the measured probability density function (e.g., as the 
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Fig. 6. Kinetic energy distribution of electrons emitted from air-exposed 
gold surfaces under UV illumination with 253.6 nm wavelength at ~ 300 K. 




Fig. 7. Illustration of an initial velocity vector vq of an emitted electron 
w.r.t. the surface normal n. 



one shown in Figure [6]), which is assigned to the corresponding 
sensor surfaces. The initial electron velocity is given as 



1^0 1 



Ekin(-^l) 



(3) 



where rrie is the electron mass at zero velocity. 

The release directions of photo-electrons follow a cosine- 
distribution |[35lll29ll . Figure [T] illustrates the initial velocity 
vector vo of an electron with respect to the surface normal 
vector n. The direction is defined by the angles and cp 
according to: 



arccos ( y 1 — 



(4) 



where X2 and X3 are uniformly distributed pseudo random 
numbers between and 1. The velocity vectors are distributed 
on a half sphere with the emitted electron in the center. 

The force acting on an electron in a time-variant electric 
field is given as 

F{r,t) = e-E{r,t) = me-r{r,t). (5) 

In Eq. [5] e is the electron charge and E (r, t) the strength 
of the electric field at electron position r and time t. The 
acceleration r (r, t) on the electron can be written as a second- 
order differential equation: 



r(r,t) 



E{f,t). 



(6) 



The electron trajectory is the solution of Eq. [6] Due to the 
inhomogeneous electric field and the complex geometry of 
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the inertial sensor, it is usually not possible to calculate an 
analytical solution. Therefore, Eq. [6] is solved numerically, 
using the explicit 4th/5th order Runge-Kutta one-step solver. 
The electric field at time t can be assumed static during 
the computation of individual electron trajectories, since their 
transit time in typical (millimeter- sized) sensor gaps is a few 
nanoseconds and the periods of the applied sinusoidal voltages 
(see Appendix |A|) are much larger. 

The local strength E {f^t) of the time- variant electric field 
inside the inertial sensor is computed using a commercially 
available finite-element software, where the same CAD model 
as used for the ray tracing has been imported for consistent 
geometry definition. It would require many calculations to 
obtain the electric fields inside the sensor for each possible 
combination of time-variant electrode voltages and test mass 
potentials that might occur during operation. Thus, the follow- 
ing approach is adopted: 

The influence of each individual body (electrodes, housing 
structure, test mass) on the total electric field is computed sep- 
arately by the finite-element software, for a constant potential 
of 1 V, while all other bodies are set to zero potential. With 
the obtained sets of electric fields, the total electric field can 
be reconstructed during electron tracing at each time step and 
for every possible combination of electrode potentials. This is 
achieved by superposition of the individual 1 V-field solutions, 
where each field is scaled according to the currently applied 
electrode voltages and the resulting test mass potential. 

The finite-element software provides the 1 V-field solutions 
as locally discretized vector fields (the meshing in the finite- 
element calculation has been selected to limit the energy 
error of each field solution below 0.1%). Therefore, during 
the numerical integration of Eq. [6j each 1 V-field has to be 
interpolated to obtain the field strength at the current electron 
position. This is done by means of a vectorial interpolation 
method based on spatial Delaunay segmentation f37| and 
bary centric-rating. Hence, the method allows to trace electrons 
at arbitrary instants of time, covering all possible electric field 
situations within the sensor. 

The photoemission of an arbitrary surface is modeled 
as the product of the quantum yield from bulk emission ^ 
(defined as the number of emitted electrons per absorbed 
photons) and the quantum yield gain which depends on 
the angle of light incidence. The gain models the influence of 
the vectorial photo-effect which might be caused by surface 
effects. The value of the gain is normalized to 1 at normal 
incidence and usually grows for increasing incidence angles 
Il29il. 

2) Electron tracer outputs: The main output are electron 
transition ratios between arbitrary sensor surfaces. The tran- 
sition ratio fi^j{t) from surface i to another surface j is 
computed as the weighted sum over all traced electrons that 
travel from surface i to surface j, divided by the weighted sum 
over all released electrons from surface i. The weight of each 
electron is given by the product (,m ' 9m ' Pm, where m is the 
unique identification number of an electron. In this context, 
is the (bulk emission) quantum yield of the emitting surface, 
Qm the quantum yield gain, and pm is the fraction of absorbed 




Fig. 8. Visualization of some electron trajectories. The trajectories are 
simulated for the LISA Pathfinder discharge system, where the test mass is 
illuminated (as shown in Figure [3} and a typical high frequency voltage for 
electrostatic sensing is present (see Eq. [Tsjin Appendix [A). The trajectories 
are valid only for a snapshot in time and notably vary due to the time-variant 
electric field. 



light at the ray impact (and electron release) position. 

Due to the applied AC voltages, the transition ratios are 
functions of time. Usually, their values significantly change 
with the amplitude and orientation of the electric field. The- 
oretically, the value at a specific point in time tx is simply 
obtained by considering only electrons that have been emitted 
exactly at time t^. Practically, the transition ratios at time tx 
are obtained by considering all simulated electrons in the small 
interval At = [tx—r < tx < tx-\-r]. The computed trajectories 
for a certain time interval At can be visualized, as shown in 
Figure [8] 

C. Analytical Discharge Model 

The simulated electron trajectories can be used to set up an 
analytical model for the instantaneous charge rate Qi^j{At) 
from an arbitrary surface i to another surface j: 



Q,^j{At) 




{At) -IvyiAt). (7) 



In Eq. [Tj M is the total number of traced electrons within 
the small time interval At and /uv(At) is the injected radiant 
UV light power. The function Si^j {At) becomes one when the 
electron m makes a transition from surface i to j within At (as 
can be identified from its computed trajectory); otherwise it is 
zero. The transition ratio fi^j {At) from surface i to surface 
j at time At is given as: 



(At) 



X^rn^l ' 9m ' Pm 



(8) 



where 5i{At) becomes one when electron m is released from 
surface i in the interval At. The instantaneous net charge rate 
of an arbitrary surface i is given as: 

N 

Q,{At) = ^ [Q^^J{At) - Qj^,{At)y (9) 
is the total number of modeled inertial sensor surfaces. 
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In order to quantify the robustness of a discharge system, 
not only the instantaneous discharge rates are of interest, but 
also the averaged discharge rates within an arbitrary time 
span T > At. The averaged electron flow from surface i to 
surface j is obtained when Eq. [7] is evaluated for all electrons 
generated within T. Moreover, the overall averaged electron 
flow Qtm^eh(^) from the TM to the EH is given by the sum 
of all individual electron flows from the test mass surfaces 
(Vz G TM) to the EH surfaces (V j G EH). Each electron flow 
from TM surface i to EH surface j is computed according 
to Eq. |7] for all electrons generated within T. The averaged 
electron flow Qeh^tm(^) from the EH to the TM is derived 
similarly. 

The individual (averaged) electron flows between test mass 
and electrode housing surfaces can be used to quantify the 
robustness AQ of a discharge system design to obtain bipolar 
test mass discharge rates (e.g., despite the presence of quantum 
yield variations between illuminated surfaces). For example, a 
positive discharge rate is obtained when: 



positive discharge rates; the feedthroughs pointing towards the 



AQ 



+ _ 



Qj 



EH^TM 



(T) 



> 1. 



(10) 



The actual value of the ratio AQ^ quantifies the robustness to 
achieve positive discharge rates; a value AQ^ 1 indicates 
that a positive discharge rate can be safely achieved, although 
the quantum yields of the EH surfaces are larger than that of 
the TM surfaces. Negative discharge rates are obtained when 



the reciprocal value of Eq. [T0| is larger than 1. The desired 
robustness gain for a negative discharge rate is AQ" ^ 1; 
hence, it can be realized even when the TM quantum yields 
are larger than that of the EH. 

IV. Analysis of the LISA Pathfinder UV Light 
Discharge System 



In this section, the discharge toolbox presented in Section III 
is used to analyze the baseline LISA Pathfinder discharge 
system 1211 122 1 as well as of a design modification that has 
been proposed to increase the system robustness. 

A. Baseline LISA Pathfinder Discharge System 

The UV light is generated by Hg-discharge lamps of which 
wavelengths shorter than 230 nm are removed by a filter so that 
only light of 253.6 nm wavelength (4.89 eV) is transmitted. 
The radiant power can be controlled by setting one of 256 
different commands, at a maximum rate of 1 Hz. As detailed 
in Appendix |A| the frequencies of the sinusoidal voltages 
applied to the different electrodes range between 60 Hz and 
100 kHz. Hence, the LPF discharge principle is classified 
"unsynchronized DC charge control", according to the def- 



inition in Section |II-C[ Measured calibration tables, relating 
the command setting and the radiant power Ijjy emitted at the 
tip of the UV light feedthroughs, are available. The dynamical 
range (defined as the ratio of maximum and minimum stable 
photon output) is approximately 100. 

The feedthrough pointing towards the TM (denoted UVLF 1 
in Figure 15 of Appendix |A|) is normally used to generate 



housing structure (denoted UVLF 2 and 3 in Figure [T5]) are 
used to generate negative discharge rates. 

The measured light distribution emitted from a feedthrough 
tip has been processed in order to obtain a two-dimensional 
distribution function (similar to the one shown in Figure [4]). 
The processed light intensity distribution is used as input 
parameter for the light source model of the developed ray 
tracing tool. The LPF inertial sensor geometry is directly 
obtained from the CAD model of the flight sensor and the 
materials with their corresponding reflection and scattering 
properties are assigned to the individual sensor surfaces. For 
the UV wavelength of 253.6 nm, the critical surface roughness 
to assume specular reflection is A/20 « 13 nm. Therefore, 
non-specular scattering effects have to be considerecj^ when 
the mean surface roughness 13 nm (which is the case 

for the electrode housing structure and the test mass comer 
spheres). 

The initial velocities of the photo-electrons are obtained 
from the kinetic energy distribution measurements described 



in Section |III-B1[ In particular, for the various gold coated 
sensor surfaces, the distribution shown in Figure [6] is used as 
input parameterization for the electron tracing tool. 

From the ray tracing and subsequent electron tracing simu- 
lations, it became clear that a large amount of photo-electrons 
do not contribute to the discharge rates as desired. Especially 
for the negative discharge rates, obtained by illumination 
of the electrode housing, only few of the electrons emitted 
inside the strongly illuminated caging finger holes reach the 
test mass. Moreover, a considerable part of the UV light is 
reflected from the EH to the TM (ptm = 11.6%), where 
72% of the emitted electrons make the transition to the 
housing (/tm^eh(^) = 0.72). The calculated transition ratios 
/eh^tm(^) and /tm^eh(^) for EH illumination with UVLF 2 
are shown in Figure |9] over one period of the injection 
voltage. The effect of the low-frequency actuation voltages 
is considered, since the shown values at each point are mean 
values which have been calculated from 10^ injection voltage 
periods within one second. 

The predicted robustness gains AQ of the baseline LPF 
discharge system to achieve bipolar test mass discharge rates 
are reported in Table |l| Assuming the quantum yields of all 
not gold coated surfaces to be zero (caging fingers, TM comer 
spheres), AQ directly shows the robustness against quantum 
yield imbalances of the EH and TM gold coatings. Thus, 
to obtain negative discharge rates, the quantum yield of the 
TM must not exceed the EH yield by the factors reported 
for AQ~ . Positive discharge rates are obtained only when the 
EH yield does not exceed the TM yield by the factors reported 
for AQ^. The robustness values in Table |l| are reported for 
three different scenarios of applied voltages: (1) no DC bias 
voltages, (2) maximum possible DC bias voltages to support 
positive test mass discharge rates, (3) maximum DC bias 
voltages to support negative discharge rates. The effect of 

^The presented analysis has been performed without characterization of 
rough sensor surfaces through BRDF measurements; however, BRDF mea- 
surements are currently being performed at PTB Braunschweig and will be 
used to refine the simulations. 
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Fig. 9. Time-variant and averaged electron transition ratios for negative 
TM discharge rates, plotted over 10 /is (corresponding to one period of the 
injection voltage). The EH and TM illumination ratios obtained with UVLF 2 
are: peh = 88.4% and ptm = 11.6% (baseline design); per = 98.4% and 
PTM = 1.6% (design modification). 



the injection and actuation voltages on the averaged transition 
ratios are considered in each scenario. 

The quantum yields of representative, air-contaminated gold 
coated samples have been observed to vary by approximately a 
factor of 3 when the sample surfaces are controlled according 
to a flight-model AIT process (uncontrolled surface samples 
have been found to vary by a factor of 70). Therefore, as 
evident from Table |lj the baseline design does not toler- 
ate the expected yield variations without application of DC 
bias voltages, even when the surfaces have been controlled 
(AQ- = 1.7). 

TABLE I 

Predicted robustness of the LISA Pathfinder discharge system 



DC Voltage Scenario 


LPF Baseline Design 


LPF Modified Design 




AQ- 




AQ- 


(1) No DC volt. 


4.3 


1.7 


4.8 


16.1 


(2) DC volt. Qtm > 


9.1 




10.6 




(3) DC volt. Qtm < 




9.8 




50.9 



B. LISA Pathfinder Discharge System Modification 

After failures to produce negative discharge rates have 
occurred during system level testing using an inertial sensor 
replica in a torsion pendulum configuration at University of 
Trento |38|, the authors investigated means to increase the 
robustness of the baseline LISA Pathfinder discharge system 
design. The developed discharge toolbox has been used the 
analyze the robustness of acceptable design modifications. The 
finally proposed modifications are summarized below: 

• Negative discharge rates: The robustness has been im- 
proved through a modification of the light injection into 
the inertial sensor. For the given constraints to leave the 
feedthrough position and the UV light source unchanged 
(inertial sensor flight hardware and mechanical interfaces 
with other parts of the LTP structure are already built), 
the optimal solution was to re-direct the light into the x- 
gap on the — z-face of the EH, such that significantly less 



light is reflected onto the TM while slightly increasing the 
electron transition ratio /eh^tm(^). The re-direction can 
be obtained by attaching a micro optical element to the 
tip of the UVLFs no. 2 and 3. 
• Positive discharge rate: Increase the gold coating of the 
strongly illuminated spherical caging finger interfaces at 
the test mass corners (see Figure [5]). The spheres have 
been partly gold coated such that the fingers do not 
damage the coating but more electrons will be released 
from the test mass (gold work function is lower than that 
of the uncoated AuPt bulk material; thus, more electrons 
will be released from the test mass). 
The optimization of the positive and negative discharge rates 
took ray-tracing and electron tracing simulations into account. 
For the negative discharge rate, the EH illumination has been 
increased from 88.4% to 98.4%. Hence, the undesired TM 
illumination has been reduced to 1.6 % (86 % less than for the 
baseline design). The transition ratios for the modification are 
also shown in Figure [9] All time-variant voltages within one 
second have been considered in the same way as described in 



Section [rV-A[ The predicted robustness values of the modified 
LPF discharge system are reported in Table [l| The achieved 
robustness values satisfy the observed quantum yield variations 
by a factor of three, also for negative discharge rates without 
application of DC bias voltages (AQ" has been increased 
from 1.7 to 16.1). 

The achieved robustness gains are still not ideal, because 
outgassing in space might produce additional surface contam- 
inations and therefore further increase the quantum yield im- 
balances. This may cause a malfunction of bipolar discharging, 
in particular when the application of DC bias voltages is not 
allowed, like for the test of continuous discharge control — 
the standard discharge mode during science measurements in 
LISA. 

V. Design and Analysis of a More Robust 
Discharge System 

Motivated by the marginal robustness of the LPF discharge 
system design, a more robust design is presented below. It can 
be directly applied to the reference inertial sensor with cubical 
test masses (see Appendix [A|. Therefore, it is proposed to be 
used for the future LISA/NGO mission, and, for any other 
space mission making use of similar inertial sensors. 

A. Performance Requirements 

In order to properly satisfy the LISA/NGO closed-loop 
charge control requirements for the fast discharge operation 
after a TM is released from the caging mechanism, and 
for continuous discharge operation during the science mea- 
surements, the discharge system shall be able to realize the 
following performance requirements. 

Assuming the expected initial test mass charge level Qtm 
according to Section |I-A[ a reasonable justification for the 
maximum required bipolar discharge rate is to discharge a 
test mass potential of ±1V in one hour. With the total test 
mass to housing capacitance of 34.2 pF 1391 , a maximum 
discharge rate requirement of ±6 • lO^e/s is obtained. A 
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reasonable minimum discharge rate to fulfill the closed-loop 
charge control requirements is ±10e/s. The required dynamic 
range (ratio between maximum and minimum discharge rate) 
for both positive and negative discharge rates is therefore 
6 • 10^ e/s. Furthermore, any value between the minimum and 
maximum discharge rate shall be realized in steps of ±10e/s 
(i.e., the charge rate resolution between the maximum positive 
and negative discharge rates corresponds to 1.2 • 10^ steps). 

B. Design of the Robust Discharge System 

From the comprehensive discharge modeling and analysis 
it becomes clear that the UV light injection of the LISA 
Pathfinder discharge system is not optimal because various 
surfaces and complex geometrical features are illuminated. 
This includes the illumination of finger holes (where the 
released electrons have poor transition ratios to reach the test 
mass), test mass corners (spherical caging finger interfaces 
are not gold coated and have no well-defined surface structure 
after de-caging), different electrodes (various disadvantageous 
electric fields have to be considered), and large parts of the 
housing structure (which can not be biased directly through 
application of DC voltages). As becomes clear from Table |lj 
the presented modifications of the LPF discharge system can 
only partially mitigate these disadvantages. 

The LISA Pathfinder UV light source hardware unit imposes 
further constraints: Only unsynchronized DC operation is 
possible and the light output performance is characterized by 
strong temperature dependency and limited lifetime of the Hg 
discharge lamps. 

From this experience, a more robust discharge system has 
been designed with the following goals: 

• Obtain sufficient robustness against quantum yield and 
work function variations such that bipolar test mass 
discharge rates can be safely achieved. 

• Satisfy the discharge rate performance requirements of 



Section |V-A| for any realistic imbalance of the surface 
photoemission properties. 
• Minimize charge noise due to discharge actuation 
The upper goals can be achieved by two major design 
changes: 

1) Change of the UV light injection into the inertial sensor 
such that the light is injected through the existing holes 
in the center of the injection electrodes which are 
mounted on the ?/-faces of the electrode housing (see 
Figure [TO]). The light shall be injected such that it is 
mainly restricted on the surfaces between the injection 
electrodes and its adjacent side on the test mass. 

2) Change the UV light source and the actuation strategy 
such that the light is only switched on when the present 
electric fields entirely suppresses the unwanted electron 
flow. 

As indicated in Section 



II-B higher UV light energies 



mitigate the effects of small work function changes due 
to molecular adsorption as well as effects related to 
surface emission. Thus, the light source characteristics of 
commercially available UV LEDs with a spectral peak 
wavelength of 240 nm (5.17 eV), a maximum output power 




Fig. 10. Proposed light injection via an UVLF, located in the existing hole of 
the injection electrode (red), on the — y-side of the electrode housing (yellow). 
Some simulated UV light rays are shown. The UVLA tip with the fiber inside 
is chamfered by 30°, such that the central light beam is deflected by 18.8°. 
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Fig. 11. Spectrum of the 240 nm UV LED as parameterized in the 
light source model of the ray tracer to generate random rays with different 
wavelengths. 



of 100 /iW, and a maximum modulation frequency of up 
to 100 MHz are assumed in the following calculations. 
The spectrum as parameterized in the light source model 
of the ray tracer is shown in Figure [TT] The effect of 
reduced wavelengths on photo-electron emission from air- 
contaminated gold surfaces is currently studied in the scope 
of an ESA technology development program 1401 . 



1) Light injection: It is proposed to inject the UV light 
through each of the two existing holes on the -^y and —y 
side of the electrode housing (see Appendix |A]). As and 
example. Figure [TO] shows the UVLF located on the —y- 



side of the electrode housing. The light injection has been 
optimized by running multiple ray-tracing simulations with 
the optimization criterion to maximize the illumination of the 
injection electrode on the ?/-face and its adjacent side on the 
test mass, while minimizing the illumination of other parts 
(e.g., the UVLF and the actuation/sensing electrodes). 

The core of the UVLF is a multimode silica fiber of 
1mm diameter with a refractive index n = 1.5 at the peak 
wavelength (same fiber as used for LPF). The fiber tip is 
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chamfered by 30° such that the central beam of the exiting 
UV Hght is deflected away from the central axis of the 
UVLF by 18.8°. The deflection is caused by the change 
of refractive index at the transition between the silica fiber 
and vacuum, as described by Snell's law. In principle, this 
method allows deflection angles up to ^ 40°, unless the 
central ray of the light beam is totally reflected inside the 
fiber. However, for the expected beam spread, parts of the 
outer beam wings are already reflected when the chamfer angle 
exceeds 30°. Refractive light deflection is used, because it 
allows to sufficiently fulfill the optimization criterion and has 
a reduced manufacturing complexity compared to a micro- 
optical element. 

From the visualization of the ray tracing simulation in 



Figure 10 it becomes clear that the light initially hits the test 
mass and is then reflected back to the injection electrode. 
Consequently, 71.5% of the light will be absorbed by the 
test mass and 23.2% by the adjacent injection electrode; 
the remaining 5.3 % are absorbed by other surfaces (mainly 
by the UVLF and the hole in the injection electrode). Note 
that the light distribution within the sensor is constrained on 
regions with simple geometry and only one dominant electric 
field (caused by the voltage applied to the injection electrode 
and by the test mass potential). 

2) Actuation strategy: Due to the obtained light distribu- 
tion, electrons are always emitted from the test mass and the 
electrode housing parts (similar to the unsynchronized DC 
charge control principle of LISA Pathfinder). Depending on 
the kinetic energy of the photo-electrons and the electric fields 
in the gap between injection electrode and test mass, they 
either reach the adjacent side or are attracted back to the 
emitting side. A robust positive test mass discharge rate is 
obtained by switching on the light only when the electric field 
is such that all electrons emitted from the EH are suppressed 
(and the transition of electrons emitted from the test mass 
is enhanced), and, vice versa for negative discharge rates. 
Figure [12] illustrates the "UV light pulsing" principle for a 
positive test mass discharge rate. The schematic shows the 
UV light injection via the chamfered UVLF, the instantaneous 
electrical field situation between the illuminated test mass and 
electrode housing surfaces, and the desired electron trajectories 
to achieve positive discharge rates. 

Since no other electrodes than the ^/-injection electrodes are 
illuminated, the already present 100 kHz sinusoidal injection 
voltage (see Eq. [13] in Appendix |A]) is adequate to properly 
control the transition of the emitted electrons from both sides. 
This requires the light pulses to be synchronized with the 
100 kHz injection voltage. The concept of the proposed UV 
light pulsing relative to a digital 100 kHz synchronization 
signal is shown in Figure [13] The digital synchronization 



signal is an input to the light source electronics and should be 
provided by the same actuation and sensing electronics which 
also generates the injection voltage. 

Furthermore, in order to describe the synchronized on/off 
switching of the UV light in a compact way, a command set 
consisting of the four parameters (n, p, o, w) at a. command 
rate of IHz is proposed for the UV light source unit (see 




Fig. 12. Schematic illustration of the pulsed-light-source design, shown for 
a positive test mass discharge rate. The light pulses are synchronized with the 
injection voltage such that emitted electrons from the EH are suppressed by 
the electric field and those from the TM are accelerated towards the EH. 
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Fig. 13. Concept of proposed UV light switching (synchronized AC charge 
control). The four actuation commands (n, p, o, w) define the number of light 
pulses per 1 Hz command period, the radiant power of a single pulse, phase 
offset of the light pulses w.r.t. the 100 kHz synchronization signal, and the 
pulse width of the light pulses. 



Figure With this concept, the load of the spacecraft 
communication bus is minimized, while achieving a high 
dynamic range of the UV output. 

The parameter n defines the number of light pulses per 1 Hz 
command period. The maximum number of light pulses is 
determined by the number of periods of the 100 kHz injection 
voltage within one second. Thus, between and 10^ light 
pulses can be generated for both the positive and the negative 
half-periods of the injection voltage. The parameter p specifies 
the radiant power of the individual light pulses. The parameter 
o sets the phase offset between the light pulses and the 100 kHz 
digital synchronization signal. Therefore, the offset parameter 
also defines whether a positive or a negative discharge rate is 
obtained (by "shifting" the light pulses to either the positive 
or the negative half-period of the injection voltage). Note that 



each of the two feedthroughs can be used to generate positive 
and negative discharge rates, providing double redundancy for 
each discharge direction. The parameter w sets the pulse width 
of the individual light pulses. 

The resolution of p, o, and w is proposed to be 8 bit, 
such that 256 steps can be commanded for each of these 
parameters. Thus, the phase offset o ranges between (i.e., 
no offset) and 255 (10 /iS offset), with a linear step size of 
^ 40 ns. By specifying the maximum pulse width w to be 
one half of the injection voltage period, the maximum pulse 
width command is 5 /iS, and the minimum possible pulse width 
command approximately 20 ns (which also defines the pulse 
width command resolution). 

The resolution of the parameters n, p, and w define the 
theoretical dynamic range of the UV light source according to 
10^ -2^ -2^ ^ 6.6-10^. This is more than 7 orders of magnitude 
larger than the dynamic range of the LISA Pathfinder UV light 
source. 

In order to obtain bipolar discharge rates, independently 
from quantum yield imbalances between EH and TM surfaces, 
the offset and the pulse duration within one half-period of 
the injection voltage are selected such that its amplitude is 
always larger than the kinetic energy of the emitted electrons. 
Assuming the maximum electron kinetic energies to be 3eV 
(which is a worst case assumption since it would correspond 
to work functions of only 2.1 eV), the phase offset and the 
pulse duration are set such that the light is only switched on 
when the voltage between the injection electrode and the TM is 
larger than 3 V. This reduces the theoretical dynamic range by 
approximately a factor of 2, but causes the undesired electron 
flow to be almost entirely suppressed. 

C. Analysis of the Robust Discharge System 

Figure [14] shows the calculated transition ratios between the 
test mass and the electrode housing surfaces over the period of 
the 100 kHz injection voltage Vinj. The effect of the actuation 
voltages is considered since the shown transition ratios are 
averaged over all injection voltage periods within one second. 
The highlighted area in the center of the first half-period 
of the 100 kHz injection voltage shows that the transition 
of electrons from the EH to the TM is almost completely 
suppressed (/eh^tm(^) ^ 0); the transition ratio from TM to 
EH is /tm^eh(^) ^ 1- Thus, the UV light is switched on 
during this time slot to obtain positive discharge rates. On the 
other hand, the electrons emitted from the test mass are almost 
completely suppressed in the center of the second half-period 
(/tm^eh(^) ^ 0); thus, the light is switched on during this 
time slot to obtain negative discharge rates. 

In both cases, the basic principle is to suppress the "un- 
wanted" electron flow as much as possible. The quantum yield 
imbalances between EH and TM surfaces merely contribute 
to the uncertainty in the amplitude of the obtained discharge 
rates, but do not influence the functionality to obtain bipolar 
discharge rates. The desired discharge rate amplitudes can be 
recovered through the amount of injected UV light power 
(assuming the dynamic range of the light source is sufficiently 
large). 



0.6 



c 0.4 

H 

0.2 




^M— >EH 
^H— >TM 



Qtm>0 



Qtm<0 



x-x^ 



0.2 



0.4 0.6 
Time [s] 



X 10 



Fig. 14. Transition ratios of the emitted electrons between EH and TM 
surfaces over the period of the injection voltage. The ideal time slots in which 
the UV light has to be switched on to obtain positive and negative test mass 
discharge rates are highlighted. 



Table |ll| shows the robustness gains of the proposed design, 
without application of DC bias voltages. When the typical 
actuation voltages for zero force/torque command (see Ap- 
pendix |A]) are considered in the analysis, bipolar discharge 
rates can be realized for quantum yield imbalances up to 
a factor of ^ 884. The effect of the actuation voltages is 
indicated by the calculated robustness gains when no actuation 
voltages are applied. 

TABLE II 

Predicted robustness of the proposed discharge system design 
No Actuation Voltages With Actuation Voltages 



DC Voltage Scenario AQ+ 



AQ- 



No DC volt. 



1238.6 



785.2 



1203.5 



Moreover, electron kinetic energies of about 3eV (corre- 
sponding to a work function of 2.14 eV) can be tolerated 
without degradation of discharge performance. This assumes 
that the phase offset o and the pulse duration w are such that 
the UV light is only switched on when the undesired electron 



flow is sufficiently suppressed (see Figure 14). 

Assuming a quantum yield imbalance of 884 (corresponding 
to smallest acceptable value of the robustness gains in Table 
when actuation voltages are applied), the discharge perfor- 
mance requirements of Section |V-A| can be fulfilled when the 
dynamic range of the light source is Iuv,max/Iuv,min = 1.6-10^. 
This range can be achieved with the presented design. 

VI. Conclusion 

A comprehensive modeling toolbox for contact-free UV 
light discharge systems has been developed and discussed. The 
results from dedicated measurement campaigns have been used 
as model inputs (e.g., width and shape of the kinetic energy 
distribution of the photo-electrons, quantum yield, reflection 
curves) to have a representative description of the underlying 
physics. The toolbox has been used to analyze the existing 
flight hardware of the LISA Pathfinder discharge system. Im- 
portant findings are that the robustness of the baseline design 
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is marginal, since it does not guarantee the capability to safely 
produce bipolar discharge rates for the observed variations 
of the measured surface properties. These findings motivated 
the definition of stringent AIT requirements, the optimization 
of the UV light injection and of the test mass coating, with 
the goal to increase the system robustness. The optimization 
has been performed by means of extensive simulations using 
the developed discharge toolbox with inputs from dedicated 
measurement campaigns. 

An analytical discharge model has been derived on the 
basis of the calculated surface illumination ratios and the 
electron trajectories. In addition to quantify the robustness of 
the LISA Pathfinder discharge system, such a model is used 
for the development of the charge control on-board software 
algorithms. Moreover, the model is the core of the discharge 
performance part of the overall LISA Pathfinder performance 
budget. 

The discharge toolbox has been used to design and optimize 
a more robust UV light discharge system, which is directly 
applicable to cubical test mass inertial sensors as used in LPF 
and planned to be used in LISA/NGO. The analysis of the 
proposed design shows a significant increase of the robustness 
and performance figures. The concept is the baseline discharge 
system design for LISA/NGO as defined in the scope of the 
LISA Mission Formulation Study BTIl . Furthermore, from the 
new design and its detailed analysis, performance requirements 
for the needed UV light source have been derived. These 
requirements are inputs to a UV light source technology 
development program, currently running under ESA contract 

ioi. 

Appendix A 
Inertial Sensor Description 

In the following, the main components of space inertial 
sensors with cubical test masses are introduced, as far as they 
are needed for the understanding and modeling of UV light 
discharge systems. 



Figure 15 shows an exploded view of the inertial sensor 
which has been used as the reference scenario for the presented 
article. Two of these sensors are part of the scientific pay load 
(LTP) on board the LISA Pathfinder spacecraft. The cubical, 
1.96 kg and 46 x 46 mm sized gold-platinum (Au73Pt27) test 
mass is nominally centered in the housing structure that carries 
the actuation/sensing and injection electrodes. All electrodes 
are surrounded by grounded guard rings to reduce fringe field 
effects. The inner parts of the sensor (except the spherical 
interfaces at the test mass corners and the caging fingers) 
are coated with a gold layer of 800 nm thickness on top of a 
200 nm thick titanium layer. Because of the large gaps between 
test mass and electrodes (4 mm along the sensitive x-axis), 
the test masses have to be caged during launch and will be 
released into free flight as soon as the spacecraft has reached 
the operational orbit. The release procedure is performed in 
two steps: First, the 8 iridium caging fingers (which apply the 
high loads during launch) are retracted along the z-axis into 
the finger holes, such that the test mass is held only by the 
two grabbing plungers. In a second step, the plungers are then 
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Fig. 15. Exploded view of the LISA Pathfinder inertial sensor including the 
cubical test mass, the electrode housing structure with the actuation/sensing, 
and injection electrodes, the caging fingers from the launch lock device, the 
grabbing plungers from the release mechanism, and the injection feedtroughs 
for the UV light optical fibers. 



simultaneously retracted such that the test mass is separated 
from their release tip^ 

On each side of the electrode housing (EH), two actua- 
tion/sensing electrodes are mounted. These electrodes are used 
to apply electrostatic forces and torques along each test mass 
degree of freedom (DoF) and to measure the linear and angular 
test mass displacement via the change of test mass to electrode 
capacitance. 

For the actuation of individual test mass DoFs, always a set 
of four electrodes is used (electrodes 1-4 for actuation of the 
test mass coordinates x and 0, 5-8 for y and 6>, and 9-12 for 
z and rj). The actuation voltages are sinusoids with different 



frequencies for different test mass DoFs. Eq. 11 shows the 
actuation voltages on the y/6>-electrodes as an example: 



Kict,5 
Kict,6 
Kict,7 
Kct,8 



-^Aiy{t) sm{27r fyt) - 
-Aiy{t) sm{27rfyt)- 
^A2y{t) cos{27r fyt) 
-A2y{t) cos{27r fyt) 



Aie{t) sm{27r fet) 
■A2e{t) cos{27r fet) 
-Aie{t) sm{27r fet) 
-A2e{t) cos{27rfet)(n) 



The voltage amplitudes Aiy and Aie are computed from the 
forces and torques as commanded by the DFACS on-board 
control software and depend on the disturbances occuring 
during spacecraft operation. The frequencies of the shown 
actuation voltages are fy = 90 Hz and fe = 180 Hz. 
The frequencies of all six actuation voltages range between 
60 — 270 Hz and are orthogonal within one DFACS actuation 

^When two metallic bodies (like the gold coated test mass and gold- 
platinum dental alloy of the release tip) have been brought into contact and 
are then separated, two opposite charges remain on the two bodies | 42 1 such 
that a random charge Qtm with arbitrary sign will be left on the test masses. 
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cycle (10 Hz), to minimize electrostatic force and torque cross- 
couplings. 

Also sinusoidal test signal voltages Vrs^i for system identi- 
fication purposes (e.g., test mass charge estimation 1281 ) and 
DC bias voltages V^ci (e.g., to assist test mass discharging) 
can be applied via the on-board software to each of the twelve 
actuation/sensing electrodes. The maximum amplitudes of the 
DC bias voltages are ±5 V. The frequency is restricted by the 
10 Hz sampling frequency of the on-board computer. In total, 
the following voltages can be applied to each actuation/sensing 
electrode: 



EH,i 



TS,i 



Vr 



DC,i 



1...12. 



(12) 



The six injection electrodes (located between the actuation 
and sensing electrodes on the y and z faces of the electrode 
housing) are used to apply high frequency sinusoidal voltages 
to bias the test mass for electrostatic sensing purposes. The 
injection voltage is equal for all six injection electrodes: 

linj,, = 4.82 • sin(2^/,nj • t) j = 1 . . . 6, (13) 

where /inj = 100 kHz. The potential of the free-floating test 
mass is given as: 

^TM = ^ (QtM + Yl C'EH,i^EH,i + ^^inj ,i Mnj ,i ) , (14) 



where Ceh,^ and Cinjj are the individual electrode capaci- 
tances, Ctot the total test mass to housing capacitance, and 
Qtm the accumulated test mass charge. 

The electric field, caused by the electrode voltages and 
by the test mass potential, affects the trajectories of the 
photo-electrons, emitted from electrode housing and test mass 
surfaces under UV illumination. The UV light is injected into 
the sensor through dedicated holes at the corners of the EH —z 
side by means of UV light feedthroughs (UVLF) with optical 
fibers inside. There are 3 UVLFs, two directed towards the 
electrode housing and one towards the test mass. The fibers 
are connected to a UV light source (fibers and light source are 



not shown in Figure 15) 
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